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ABSTRACT 



^ ' Context. There is still no consensus about progenitor masses of Type HP supemovae. 

Aims. We study a normal Type IIP SN 1999em in detail and compare it to a peculiar Type IIP SN 1987A. 

Methods. We computed the hydrodynamic and time-dependent atmosphere models interpreting simultaneously both the photometric and 
spectroscopic observations. 

• Results. The bolometric light curve of SN 1999em and the spectral evolution of its Ha line are consistent with a presupernova radius of 
\^ ' 500 + 2OO/?0, an ejecta mass of 19.0 ± 1.2M0, an explosion energy of (1.3 ± 0.1) x 10^' erg, and a radioactive '^Ni mass of 0.036 ± 0.009Mo. 
A mutual mixing of hydrogen-rich and helium-rich matter in the inner layers of the ejecta guarantees a good fit of the calculated light curve to 
that observed. Based on the hydrodynamic models in the vicinity of the optimal model, we derive the approximate relationships between the 
basic physical and observed parameters. The hydrodynamic and atmosphere models of SN 1999em are inconsistent with the short distance of 
7.85 Mpc to the host galaxy. 

O Conclusions. We find that the hydrogen recombination in the atmosphere of a normal Type IIP SN I999em, as well as most likely other Type 
^ IIP supernovae at the photospheric epoch, is essentially a time-dependent phenomenon. It is also shown that in normal Type IIP supernovae the 
CZ3 homologous expansion of the ejecta in its atmosphere takes place starting from nearly the third day after the supernova explosion. A comparison 



Oh 



of SN 1999em with SN 1987A reveals two very important results for supernova theory. First, the comparability of the helium core masses and 
the explosion energies implies a unique explosion mechanism for these core collapse supernovae. Second, the optimal model for SN 1999em is 
characterized by a weaker '^Ni mixing up to x 660 km s"* compared to a moderate '*Ni mixing up to ~ 3000 km s"' in SN 1987A, hydrogen 
being mixed deeply downward to ~ 650 km s"' . 
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1. Introduction 

The supernova (SN) 1999em was discovered by the Lick 
Observatory Supernova Search on October 29.44 UT in the 
nearly face-on SBc galaxy NGC 1637 (Li 1999). Detected 
shortly after the explosion at an unfiltered magnitude of ~ 13.5, 
SN 1999em was bright enough to be observed well both pho- 
tometrically and spectroscopically for more than 500 days 
(Hamuy et al. 2001; Leonard et al. 2002; Elmhamdi et al. 
2003). SN 1999em was the first Type Il-plateau supernova (SN 
IIP) detected at both X-ray and radio wavelengths, being the 
least radio luminous and one of the least X-ray luminous SNe 
(Pooley et al. 2002). The X-ray data indicated the interaction 
between SN ejecta and a pre-SN wind with a low mass-loss 
rate of ~ 2 x lO'^MQyr\ Leonard et al. (2001) presented the 
first spectropolarimetry of SN IIP based on the optical observa- 
tions of SN 1999em during ~ 160 days after SN discovery. The 
weak continuum polarization increasing from p « 0.2% on day 
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1 to p w 0.5% in the final observations was detected with an 
unchanging polarization angle x 160°. To date this event has 
become the most studied normal SN IIP. 

In order to study properly any individual object, an accurate 
distance is required in addition to high quality observations. It 
is this point that is not clear for SN 1999em. There was a rea- 
sonable agreement in determining the distance to SN 1999em 
using the expanding photosphere method (EPM; Kirshner & 
Kwan 1974): 7.5+0.5 Mpc (Hamuy et al. 2001), 8.2 + 0.6 Mpc 
(Leonard et al. 2002), and 7.838 ± 0.331 Mpc (Elmhamdi et al. 
2003). However, Leonard et al. (2003) identified 41 Cepheid 
variable stars in the galaxy NGC 1637 and found that the 
Cepheid distance to the host galaxy was 11.7 ± 1.0 Mpc. This 
distance to SN 1999em is ~ 50% lai-ger than the earlier EPM 
distance estimates. Assuming a correlation between the total 
^^Ni mass and the explosion energy for SNe IIP, Nadyozhin 
(2003) evaluated the distance to the galaxy NGC 1637 of 1 1.08 
Mpc close to the Cepheid distance. Using the spectral-fitting 
expanding atmosphere model (SEAM), Baron et al. (2004) ob- 
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tained the distance to SN 1999em of 12.5 + 1.8 Mpc in a good 
agreement with the Cepheid distance scale. Finally, study- 
ing various ingredients entering the original EPM, Dessart & 
HilUer (2006) improved this method and also achieved better 
agreement to the Cepheid distance with an estimate of 1 1 .5 ± 1 .0 
Mpc. 

Starting the investigation of SN 1999em, a normal SN IIP, 
it is impossible to pass by the well-known and studied SN 
1987A in the Large Magellanic Cloud (LMC), a peculiar SN 
IIP. Recent progress in reproducing the bolometric light curve 
observed in SN 1987 A with a modern hydrodynamic model 
(Utrobin 2004) and in modelUng the Ha profile and the Ba 
II 6142 A Une in SN 1987A at the photospheric phase with 
a time-dependent approach for the SN atmosphere (Utrobin & 
Chugai 2002, 2005) makes this experience very instructive for 
other SNe IIP. First, it turned out that a good agreement be- 
tween the hydrodynamic models and the photometric observa- 
tions of SN 1987A did not guarantee a correct description of 
this phenomenon as a whole. Second, the strength of the Ho; 
line and its profile provided hard constraints on the hydrody- 
namic models. The most important lesson from this study of 
SN 1987A is that we have to take both the photometric and 
spectroscopic observations into account to obtain an adequate 
hydrodynamic model (Utrobin 2005). Now such an approach 
should be appUed to SN 1999em. 

Here we present the comprehensive hydrodynamic study of 
SN 1999em complemented by the atmosphere model with the 
time-dependent kinetics and energy balance. A brief descrip- 
tion of the hydrodynamic model and the atmosphere model 
based on it is given in Sect. 2. The study of SN 1999em we 
begin with a construction of the optimal hydrodynamic model 
(Sect. 3). Then we continue with a question: at what distance 
to SN 1999em, short or long, the hydrodynamic model and the 
proper atmosphere model are consistent with the photometric 
and Ha line observations (Sect. 4). The time development of 
the optimal hydrodynamic model is presented in Sect. 5, and its 
general regularities in Sect. 6. The basic relationships between 
the physical and observed parameters for SNe IIP similar to 
the SN 1999em event are obtained in Sect. 7, while in Sect. 8 
we address the comparison of SN 1999em with SN 1987A. In 
Sect. 9 we discuss our results from the theoretical and obser- 
vational points of view. Finally in Sect. 10 we summarize the 
results obtained. 

We adopt here a distance of 11 .7 Mpc (Leonard et al. 2003), 
a recession velocity to SN 1999em of 800 kms ' (Leonard 
et al. 2002), an explosion date of ID 2451476.77, and a total 
extinction Ay = 0.31 (Baron et al. 2000; Hamuy et al. 2001; 
Elmhamdi et al. 2003). 

2. Supernova modelling and input physics 

Keeping in mind the importance of a hydrodynamic study and 
of atmosphere modelling, including the time-dependent kinet- 
ics and energy balance for the interpretation of the SN 1987A 
phenomenon, we use this approach to investigate SN 1999em. 

A hydrodynamic model is computed in terms of radiation 
hydrodynamics in the one-group approximation taking into ac- 
count non-LTE effects on the average opacities and the thermal 



emissivity, effects of nonthermal ionization, and a contribution 
of lines to opacity as in the case of SN 1987A (Utrobin 2004). 
Note that the bolometric luminosity of an SN is calculated by 
including the retardation and limb-darkening effects. 

The atmosphere model includes the time-dependent ioniza- 
tion and excitation kinetics of hydrogen and other elements, 
the time-dependent kinetics of molecular hydrogen, and the 
time-dependent energy balance (Utrobin & Chugai 2005). The 
density distribution, chemical composition, radius of the pho- 
tosphere, and effective temperature are provided by the corre- 
sponding hydrodynamic model. The obtained time-dependent 
structure of the atmosphere is then used to calculate synthetic 
spectra at selected epochs. The spectra are modelled by the 
Monte Carlo technique suggesting that the photosphere diffu- 
sively reflects the incident photons and that the line scattering is 
generally non-conservative and is described in terms of the line 
scattering albedo. The Thomson scattering on free electrons 
and Rayleigh scattering on neutral hydrogen are also taken into 
account. 

3. Optimal hydrodynamic model 

Elmhamdi et al. (2003) have constructed the "UBVRI" bolo- 
metric light curve of SN 1999em from the corresponding pho- 
tometric data. To account for a possible contribution of the 
missing infrared bands, we add a value of 0.19 dex, taken from 
Eknhamdi et al. (2003), to the "UBVRI" luminosity and adopt 
the resultant Ught curve as the bolometric Ught curve of SN 
1999em. Our aim is to find an adequate hydrodynamic model 
that reproduces the photometric and spectroscopic observations 
of SN 1999em. To fit the bolometric light curve of SN 1999em, 
various hydrodynamic models were explored. The bolometric 
light curve of SN 1999em is fitted by adjusting the pre-SN 
radius Rq, the ejecta mass Me„y, and the explosion energy E, 
along with the density distribution in the pre-SN model and 
its chemical composition in the transition region between the 
hydrogen-rich envelope and metal/heUum core. The radioactive 
^^Ni mass is reliably measured by the light curve tail after the 
plateau phase. As the pre-SN model of SN 1999em, we con- 
sider non-evolutionary models similar to those of SN 1987A 
(Utrobin 2004), but for the outer layers assume the standard so- 
lar composition with the mass fraction of hydrogen X = 0.735, 
heUum Y = 0.248, and the metallicity Z = 0.017 (Grevesse & 
Sauval 1998) taking a normal spiral nature of the host galaxy 
NGC 1637 into account. The best version of such a fitting was 
obtained with the optimal model that is characterized by the 
basic parameters: the pre-SN radius of 500 Rq, the ejecta mass 
of 19 Mq, and the explosion energy of 1.3 x 10^' erg (model 
Dll in Table 1). 

The density profile of the pre-SN model consisting of a cen- 
tral white dwarf like core and an outer envelope with the size 
of the red supergiant is shown in Fig. 1. In the calculations, 
the 1.58 Mq central core is removed from the computational 
mass grid and assumed to collapse to a neutron star, while the 
rest of the star is ejected by the SN explosion, which is mod- 
elled by a piston action near the edge of the central core. The 
pre-SN model has a heterogeneous chemical composition con- 
taining a 5.6 Mq helium core and a 11.9 Mq outer shell of 
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Fig. 1. Density distribution with respect to interior mass (a) and radius 
(b) for the pre-SN model D 1 1 . The central core of 1 .58 M© is omitted. 
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Fig. 2. The mass fraction of hydrogen {solid line), heUum (long dashed 
line), heavy elements (short dashed line), and radioactive '*Ni (dotted 
line) in the ejecta of model Dl 1. 



the solar chemical composition (Fig. 2). Note that the helium 
cores up to about 8 Mq are consistent with the observations 
(Sect. 6.2). There is no sharp boundary between the hydrogen- 
rich and helium-rich layers in the ejecta of the optimal model. 
Hydrogen-rich material is mixed into the central region, and 
helium-rich material, in turn, is mixed outwards. It is evident 
that such a distribution of hydrogen and helium implies a strong 
mixing at the helium/hydrogen composition interface. The fact 
that the radioactive ^^Ni is confined to the innermost ejected 
layers suggests its weak mixing during the SN explosion. 

In Fig. 3 we show the very good match between the bolo- 
metric light curve calculated for the optimal hydrodynamic 
model and the one observed for SN 1999em (Elmhamdi et 
al. 2003). Note that hereafter tots is the time in the observer's 
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Fig. 3. Comparison of the calculated bolometric light curve of model 
Dll (solid line) with the bolometric data of SN 1999em obtained by 
Elmhamdi et al. (2003) (open circles). 

Table 1. Hydrodynamic models for distances of 7.85 and 1 1.7 Mpc. 
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frame of reference. The model agrees well with the observed 
tail of the bolometric light curve for the total ^^Ni mass of 
0.036 Mq, the bulk of the radioactive ^''Ni being mixed in the 
velocity range < 660 km s"^ (Fig. 10) in order to reproduce the 
observed transition from the plateau to the tail. 



4. Hot profile: evidence against short distance 

It is quite clear that any well-observed SN should be described 
by a unique hydrodynamic model in combination with the at- 
mosphere model based on it. The experience during the study 
of SN 1987A showed that such a combination of the models 
had to fit not only the photometric, but also the spectroscopic 
observations. We believe that the adequate hydrodynamic and 
atmosphere models of SN 1999em are able to distinguish be- 
tween the short and long distances. A distance of 7.85 Mpc, the 
average value of the EPM distance estimates, is taken as the 
short distance, and the Cepheid distance of 11.7 Mpc is taken 
as the long distance. 

As demonstrated above, we constructed the optimal hydro- 
dynamic model for the long distance (model Dll in Table 1). 
The model reproduces the observed bolometric light curve of 
SN 1999em very well (Fig. 4a). In Fig. 4b the calculated ex- 
pansion velocity at the photosphere level, the photospheric ve- 
locity, is compared with the radial velocities at maximum ab- 
sorption of the different spectral lines measured by Hamuy et 
al. (2001) and Leonard et al. (2002). The photospheric velocity 
of model Dll is consistent with the observed points, at least 
for the first 60 days. To verify the hydrodynamic model by 
matching the constraints from the spectral observations of SN 
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Fig. 4. Hydrodynamic model Dl 1 and Ho- line for the distance of 1 1.7 
Mpc. Panel (a): the calculated bolometric light curve {solid line) com- 
pared with the observations of SN 1999em obtained by Elmhamdi et 
al. (2003) (open circles). Panel (b): calculated photospheric velocity 
{solid line) and radial velocities at maximum absorption of spectral 
lines measured by Hamuy et al. (2001) and Leonard et al. (2002) {open 
circles). Panel (c): Ha profiles, computed with the time-dependent 
approach {thick solid line) and with the steady-state model {dotted 
line), overplotted on the observed profile on day 26.24, as obtained 
by Leonard et al. (2002) {thin solid line). Panel (d): the same as panel 
(c) but for day 52.14. 



1999em, we examined the Ha profile on days 26.24 and 52.14 
computed in the time-dependent approach. 

In the time-dependent atmosphere model for SN 1 987 A, we 
considered two extreme cases to allow for the uncertainty of 
our approximation in a description of the ultraviolet radiation 
field: (i) the photospheric brightness is black-body with the ef- 
fective temperature (model A); (ii) the photospheric brightness 
corresponds to the observed spectrum (model B) (Utrobin & 
Chugai 2005). For SN 1999em the photospheric brightness in 
model B is black-body with the effective temperature and the 
corresponding brightness reduction taken from the SN 1987A 
observations. 

The time-dependent approach with model B satisfactorily 
reproduces the strength of the Ha emission component on day 
26.24 with some emission deficit near the maximum in com- 
parison to what was observed (Fig. 4c). This deficit is not sig- 
nificant because model A, calculated for this phase and not 
plotted on Fig. 4c for the sake of clarity, gives the emission 
component with a relative flux of 3.1 at the maximum, which 
is much higher than the observed one, so the real situation is 
somewhere between these two cases and is closer to model B. 
Note that in SN 1987 A the Ha emission component demon- 
strated the same behavior in the early phase (Utrobin & Chugai 
2005). In contrast, the Ha absorption component calculated in 



Fig. 5. Hydrodynamic model D07 and Ha line for the distance of 7.85 
Mpc. See the Fig. 4 legend for details. 



model B is stronger than that observed on day 26.24 in SN 
1999em. This discrepancy is presumably related to a poor de- 
scription of the ultraviolet radiation at frequencies between the 
Balmer and Lyman edges. This radiation interacts with numer- 
ous metal fines and controls the populations of hydrogen lev- 
els. It is very important that both the emission and absorption 
components of the Ha line calculated in model B extend over 
the whole range of the observed radial velocities from -15 000 
km s~^ to 15 000 km s~^ (Fig. 4c). Unfortunately, the calculated 
absorption component runs above the observed one in a radial 
velocity range between -15 000 kms ' and -12500 kms 
On day 52.14 the emission component of the Ha fine com- 
puted with the time-dependent approach with model B fits the 
observed one fairly well, while the absorption component is 
still stronger than the observed one (Fig. 4d). Thus, we may 
state that the above hydrodynamic and atmosphere models are 
in a good agreement with the photometric and spectroscopic 
observations of SN 1999em. 

Now let us pay attention to the time-dependent effects in 
hydrogen lines of SN 1999em, a normal SN IIP. The time- 
dependent approach with model B reproduces the strength of 
the Ha line at day 26.24 and day 52.14 fairly well (Figs. 4c 
and 4d). In contrast, a steady-state model B demonstrates an 
extremely weak Ha line on days 26.24 and 52.14. This reflects 
the fact that the steady-state ionization is significantly lower 
than in the time-dependent model, while the electron tempera- 
ture is too low for the coUisional excitation of hydrogen. The 
radioactive ^^Ni is mixed too weakly to affect the ionization and 
excitation of hydrogen and other elements in the atmosphere at 
the plateau phase. Thus, it is possible to conclude that the hy- 
drogen recombination in the atmosphere of SN 1999em during 
the whole plateau phase is essentially a time-dependent phe- 
nomenon. 
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Fig. 6. Evolution of a normal SN IIP illustrated by the bolometric light 
curve of the optimal model. Lettered dots mark the specific times in 
the evolution. The end time of the shock breakout phase, fi, is shown 
in Fig. 8a. The cross notes the moment of complete hydrogen recom- 
bination. The numbers indicate the total optical depth of the ejecta at 
the corresponding times. 
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In turn, the short distance results in the hydrodynamic 
model with the following parameters: the pre-SN radius of 375 
Rq, the ejecta mass of 16 Mq, and the explosion energy of 
6.86 X 10^" erg (model D07 in Table 1). In this model the cen- 
tral core of 1.4 Mq is assumed to collapse to a neutron star. 
The hydrodynamic model fits the observed bolometric light 
curve of SN 1999em but for the ^^Ni mass of 0.0162 Mq, the 
most of ^^Ni being mixed in the velocity range < 580 kms ' 
(Fig. 5a). The photospheric velocity curve runs well below the 
observed points (Fig. 5b), and it might be expected that the 
spectral lines in this model would be narrower than those ob- 
served. Indeed, on days 26.24 and 52.14 both the emission and 
absorption components of the Ha line computed with the time- 
dependent approach with model B are significantly narrower 
than those observed in SN 1999em (Figs. 5c and 5d). The hy- 
drodynamic model D07 agrees fairly well with the observed 
bolometric light curve, but the relevant atmosphere model fails 
to reproduce the Ha profile observed in SN 1999em. We thus 
conclude that the short distance of 7.85 Mpc should be dis- 
carded. 

5. Development of the optimal model 

Although the major issues of Ught curve theory were recog- 
nized a long time ago (Grassberg et al. 1971; Falk & Arnett 
1977), it is useful to consider an SN outburst in a more detailed 
approach. This study of SN 1999em provides a good oppor- 
tunity to examine the time development of a normal SN IIP. 
Figures 6 and 8a show the following stages in the observer 
time scale: a shock breakout (tabs ^ ^i). an adiabatic cooUng 
phase (ti < tabs ^ h), a phase of cooling and recombination 
wave (CRW) (/2 < tabs ^ h), a phase of radiative diffusion 
cooUng (^3 < tabs ^ U), an exhaustion of radiation energy 
(t4 < tabs ^ ts), a plateau tail (ts < tabs ^ t(,), and a radioactive 
tail (tabs > te)- In addition to the above stages there are two 
specific points: a transition from acceleration of the envelope 



Fig. 7. Propagation of the shock through the ejected envelope of model 
Dll. Velocity profiles with respect to interior mass are plotted at t = 
0.0065, 0.0558, 0.1199, 0.2416, 0.6234, 0.8937 days, and f > 10 days 
after the SN explosion. 

matter to a homologous expansion and a moment of complete 

hydrogen recombination. In the optimal model Dll, the char- 
acteristic moments are ti ^ 0.93 days, r2 ~ 18 days, ^ 94 
days, ?4 « 116 days, ts « 124 days, and ~ 150 days, and the 
complete hydrogen recombination occurs at fe ~ 1 11. 3 days. 
Below we consider the basic stages in the time development of 
the optimal model. 

5.1. Shock breakout 

The explosion of the star is assumed to be triggered by a piston 
action near the edge of the central core immediately after the 
epoch zero, t = 0. From here on, t is the time in the comov- 
ing frame of reference. This energy release generates a strong 
shock wave that propagates towards the stellar surface. In mov- 
ing out of the center, the shock wave heats matter and acceler- 
ates it to velocities increasing outward and exceeding the local 
escape velocity. From t - 0.0065 days to f = 0.2416 days, the 
shock wave, propagating outward from the compact dense core 
of the pre-SN (Fig. 1), is attenuated slightly due to the spher- 
ical divergence (Fig. 7). It then reaches the outermost layers 
with a sharp decUne in density and after t = 0.6234 days gains 
strength and accelerates due to the effect of hydrodynamic cu- 
mulation. Only a small portion of the star undergoes this accel- 
eration and acquires a high velocity. For example, the layers of 
velocities exceeding 5000 km s"' have a mass of w O.478M0. 

By day 0.8651 the shock wave reaches the stellar surface 
and then begins to heat the external layers, so that the color 
temperature jumps to 3.84 x 10^ K at day 0.8823, the effective 
temperature increases up to 1.76 x 10^ K at day 0.8943, and 
the bolometric luminosity rises accordingly up to 6.46 x lO''^ 
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Fig. 8. Shock breakout in model Dll. Arrow indicates the moment of 
0.8803 days when the velocity at the stellar surface reaches escape ve- 
locity. Panel (a): the calculated bolometric light curve and the lettered 
dot marking the end time of the shock breakout phase. Panel (b): the 
calculated photospheric radius (solid line), the effective temperature 
(dashed line), and the color temperature (dotted line) as a function of 
the observer's time. 

ergs"^ at day 0.8943 (Fig. 8). Note that the maximum of the 
color temperature coincides closely with the moment of 0.8803 
days when a velocity at the stellar surface reaches the escape 
velocity. 

A very rapid rise in the bolometric luminosity to maximum, 
starting at day 0.8651, instantly changes a growth of the pho- 
tospheric radius into its reduction because of the intense ra- 
diative losses of energy in the outermost layers (Fig. 8b). At 
day 0.8691 these layers begin to move outward, and the addi- 
tional coohng by adiabatic expansion makes the reduction of 
the photospheric radius more noticeable. At the same time the 
envelope expansion is favorable to a photon diffusion, decreas- 
ing the characteristic diffusion time. The photon diffusion even- 
tually dominates, stops this reduction at day 0.8802, and then, 
along with the envelope expansion, blows the photospheric ra- 
dius outward. 

When a velocity at the stellar surface exceeds the escape 
velocity, the outside layers of the star begin to cool rapidly, the 
color temperature begins to decrease from its maximum value, 
and both the effective temperature and the luminosity decrease 
somewhat later. A narrow peak in the bolometric luminosity 
forms as a result (Fig. 8a). The peak has a width of about 0.02 
days at a half level of the luminosity maximum. Most of its 
radiation is emitted in an ultraviolet flash. The total number 
of ionizing photons above 13.598 eV for the whole outburst is 
2.768 X 10^^. During the first 1.234 days the number of ioniz- 
ing photons is 90% of the total number, and the radiated energy 
adds up to 1.59 x 10*^ ergs. We conditionally define the transi- 



tion time between the shock breakout and the adiabatic coohng 
phase, ti , as the time at which the bolometric luminosity drops 
by one dex from its maximum value, with the understanding 
that the adiabatic cooling becomes essential soon after the on- 
set of the envelope expansion. In the optimal model this transi- 
tion time is nearly 0.93 days. 

Scattering processes are fundamentally nonlocal and re- 
sult in exceeding the color temperature over the effective one 
when they are dominant (Mihalas 1978; Sobolev 1980). During 
the shock breakout, the adiabatic cooling phase and the CRW 
phase scattering processes dominate those of true absorption 
near the photospheric level in the ejecta of the optimal model 
and, as a consequence, lead to a color temperature that is 
considerably higher than the eff"ective temperature. Figure 8b 
shows that the maximum of the color temperature coincides 
very closely in time with the local minimum of the photo- 
spheric radius at day 0.8802 and comes before the maximum 
of the effective temperature and the bolometric luminosity. 

5.2. Adiabatic cooling phase 

During the shock passage throughout the star, the gas is heated 
up to about 10^ K and, as a consequence, is totally ionized. 
Both the internal gas energy and the radiation energy increase 
in the SN envelope. After the shock breakout, a dominant pro- 
cess in the subphotospheric, optically thick layers of the ex- 
pelled envelope is the cooling by adiabatic expansion. The adi- 
abatic losses drastically reduce the stored energy and com- 
pletely determine the evolution of the bolometric luminosity 
during the adiabatic cooling phase (Figs. 6 and 8a). Such be- 
havior of the luminosity lasts till the gas and radiation temper- 
atures drop to the critical value in the subphotospheric layers, 
and a recombination of hydrogen, the most abundant element in 
the ejected envelope, starts. In these layers hydrogen becomes 
partially ionized by around day 18. For the adiabatic cooling 
phase, the color and effective temperatures drop rapidly from 
1.35 X 10^ K and 9.32 x lO'^ K, respectively, at day 0.93 to 
6650 K and 6560 K at day 18. Their ratio also reduces as the 
contribution of scattering processes to the opacity decreases. 

In an optically thick medium, the strong shock wave prop- 
agates almost adiabatically throughout the stellar matter. When 
the shock wave emerges on the pre-SN surface, the adiabatic 
regime is broken and transforms into the isothermal regime. 
This transformation takes place during the adiabatic cooling 
phase and gives rise to a thin, dense shell at the outer edge of 
the ejected envelope. The dense shell is arising from day 0.945 
to day 1.022, reaching the density contrast of ~ 210 at a veloc- 
ity of ~ 12 300 kmc The formation of the dense shell starts 
in the optically thick medium at the optical depth of ~ 8, and 
ends in semi-transparent medium at the optical depth of ~ 1 . By 
day 1 8 the shell accelerates to a velocity of ~ 1 3 400 km c"' re- 
ducing the density contrast to a value of ~ 2.5. The matter in the 
dense shell is subject to the Rayleigh-Taylor instabiUty, which 
may result in the strong mixing of matter (Falk & Arnett 1977). 
The latter, in turn, may prevent the thin shell-like structure from 
forming in the outer layers of the envelope. 
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Fig. 9. Time dependence of the velocity of different mass shells (solid 
line) in the ejected envelope of model DU and the photospheric ve- 
locity (dashed line). Full and open circles indicate moments when a 
velocity of the mass shell reaches 99% and 99.9%, respectively, of its 
terminal velocity. The arrow marks the mass shell that has undergone 
an additional acceleration due to the expansion opacity. 

5.3. Homologous expansion 

The SN explosion causes an acceleration of the envelope mat- 
ter. At late times, when the acceleration of the ejecta be- 
comes negligible, the envelope matter expands homologously. 
Evidently, there is a transition from the acceleration to the ho- 
mologous expansion. This transition for different mass shells 
occurs at different times: the deeper the layer, the later the tran- 
sition (Fig. 9). Thus, the transition time of the whole expelled 
envelope is determined by the deepest layer. If we consider 
physical processes in an SN atmosphere, the relevant transi- 
tion time is given by the photospheric location. It is clear that 
the transition time for the SN atmosphere is much shorter than 
for the whole envelope. For instance, the atmosphere layers ex- 
ceed 99% and 99.9% of their terminal velocities starting from 
nearly day 2.8 and day 11, respectively (Fig. 9). 

It is interesting that a monotonic increase in the transition 
time with embedding inward the ejecta is broken at the 99.9% 
level of the terminal velocity. There is a prominent feature at a 
velocity of ~ 14 500 km c ' marked in Fig. 9. As we show later 
in Sect. 6.5, it is the result of an additional acceleration due to 
the resonance scattering in numerous metal lines. 

Homologous expansion is characterized by a density dis- 
tribution frozen in velocity space and scaled in time as oc r^. 
Such a density profile for the optimal model is shown at f = 50 
days in Fig. 10. In the inner region of the ejected envelope there 
is a dense shell with a density contrast of ~ 20 produced by the 
^^Ni bubble at a velocity of ~ 660 kmc"^ The density dis- 
tribution above the ^^Ni bubble shell is nearly uniform up to 
~ 2000 kmc"'. The outer layers with velocities in the ranges 



Fig. 10. The density and the ^''Ni mass fraction as a function of the 
velocity for model Dll at f = 50 days. Dot-dash lines are the density 
distribution fits p oc v"' and p oc v~'^. 



3000-6900 km s^^ and 6900-13000 km s ' may be fitted by 
an effective index n = -51np/51n r of 7 and 9.5, respectively, 
as seen from the density distribution with respect to the expan- 
sion velocity (Fig. 10). At a velocity of ~ 13 000 kmc there 
is another dense shell with a density contrast of ~ 3 originated 
during the adiabatic cooling phase. Note that both dense shells 
slightly change their density profiles with time. For example, 
the evolution of the ^*Ni bubble shell is clearly demonstrated 
by Fig. 11a. 

5.4. Cooling and recombination wave 

As the envelope expands, cooling by radiation - a cooling and 
recombination wave - occurs and completely dominates the lu- 
minosity of the SN by about day 18. From this time to nearly 
day 94 the bolometric luminosity is mainly determined by 
properties of the CRW (for details see Grassberg & Nadyozhin 
1976; Imshennik & Nadyozhin 1989). A main property of the 
CRW consists in generating virtually the entire energy flux car- 
ried away by radiation within its front. This radiation flux ex- 
hausts the thermal and recombination energy of cooling matter 
by means of recombination. At the inner edge of the CRW, the 
radiation flux is negligibly small and matter is ionized, but at 
the outer edge the flux is equal to the luminosity of the star, and 
matter completely recombines. 

During the CRW phase the evolution of the density in the 
ejecta illustrates the homologous expansion (Fig. 11a). The be- 
havior of the gas and radiation temperatures reflects two impor- 
tant facts (Fig. lib). First, in the subphotospheric layers in an 
optically thick medium, under conditions very close to LTE, the 
radiation temperature is virtually equal to the gas temperature. 
Second, in the transparent layers, beyond the photosphere, the 
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Fig. 11. Evolution of the density (a) and the gas (solid line) and radi- 
ation (dotted line) temperatures (b) as a function of mass for model 
Dll. Profiles are plotted at t = 10, 20, 40, 60, 80, 100, and 120 days. 
Star indicates the position of the photosphere. 
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Fig. 12. Evolution of the interior luminosity (a) and the fraction of 
ionized hydrogen (b) as a function of mass for model Dll. Profiles 
are plotted at t = 10, 20, 40, 60, 80, 100, and 120 days. Star indicates 
the position of the photosphere. 



teristics described above are evident (Figs. 12a and 12b). Hence 
in this period the CRW occurs. 

From t = 20 days to t = 100 days, the CRW propagates 
through a mass of ^ 15.96Mq, while the outermost layers, 
those responsible for the luminosity during the first 18 days, 
contribute only about 0.1 3M0. It is interesting that for this pe- 
riod the gas and radiation temperatures at the photosphere po- 
sition are nearly constant in time and equal to about 5600 K 
(Fig. lib). And from t = 40 days to t = 100 days, the photo- 
sphere is located at the outer edge of the hydrogen recombina- 
tion front (Fig. 12b). 

After f = 100 days the radiation flux at the inner edge of the 
CRW becomes comparable to the luminosity of star (Fig. 12a). 
This effect is caused by radiative diffusion from the central en- 
velope layers where the energy of radioactive decays is being 
released. Later on the luminosity is entirely determined by ra- 
diative diffusion and not by the CRW. 

5.5. Radiative diffusion cooling 

A cooling by radiative diffusion starts in the optically thick 
expelled envelope at day 94 and ends in the semi-transparent 
medium at day 116 (Fig. 6). Shortly before the end of the phase 
of radiative diffusion cooling, the complete hydrogen recombi- 
nation occurs at ?H ~ 111.3 days. In the optimal model the ra- 
diative diffusion takes place in the inner layers with the nearly 
uniform density distribution and results in the characteristic 
shoulder of the bolometric luminosity. It should be noted that 
this behavior of the light curve was discussed by Falk & Arnett 
(1977). To describe photon diffusion, we may use a simple one- 
zone approximation (Amett 1979), ignoring the contribution of 
radioactive decays to the bolometric light curve: 

d In Lboi 1 nRd 



dt 



tdif 



tdif 



3 c 



■Td , 



(1) 



where tdif is the diffusion time, Rd the characteristic radius of 
the diffusion region, and Td the optical depth of this region. 

Evidently, it is possible to estimate the diffusion time from 
both the decline of the bolometric light curve and the physical 
properties of the difliision region. The decline of the luminos- 
ity shoulder between day 94 and day 116 measures a diffusion 
time of ^ 35 days (Fig. 6). On the other hand, the radius of the 
diffusion region of « 1 .2 x 10^^ cm and its optical depth of » 63 
evaluated at day 105 give a diffusion time of x 30 days. These 
estimates are consistent with the above description of photon 
diffusion (1). Thus, the luminosity shoulder at the phase of ra- 
diative diffusion cooling is approximated well by the simple 
model with the constant diffusion time. Note that a transition 
from the constant diffusion time to the variable diffusion time, 
which is getting shorter, transforms the linear shoulder of the 
luminosity in logarithmic scale into a convex light curve at the 
end of the plateau. 



gas temperature turns out to be lower than the radiation tem- 
perature because of the weak interaction between matter and 
radiation and of the dominant role of adiabatic losses. From an 
inspection of the interior luminosity and the fraction of ionized 
hydrogen from t = 20 days to ^ = 100 days, the CRW charac- 



5.6. Exhaustion of radiation energy 

During the entire outburst, the total radiation energy in the 
ejecta is much greater than its total internal gas energy. In the 
inner optically thick layers of the envelope, in which radiation 
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Fig. 13. The plateau tail of the bolometric light curve in the optimal 
model. Panel (a): calculated bolometric light curve (thick solid line) is 
compared with the gamma-ray luminosity (thin solid line). Panel (b): 
evolution of the interior luminosity as a function of velocity. Profiles 
are plotted at t = 126, 130, 134, 142, 146, 150, 154, and 162 days. 



Table 2. Hydrodynamic models based on the optimal model Dl 1. 



Model 


Remarks 


Dec 


no dense central core (Fig. 14a) 


Ldn 


less dense outer layers (Fig. 14c) 


Mdn 


more dense outer layers (Fig. 14c) 


Hoi 


helium-rich composition in the outer layers: 




X = 0.368, Y = 0.615, andZ = 0.017 (Fig. 15a) 


Zol 


low metallicity in the outer layers: 




X = 0.735, Y = 0.259, and Z = 0.006 (Fig. 15c) 


Sci 


sharp mctals/He/H composition interface (Fig. 16a) 


Hec 


helium core of about 8 M© (Fig. 16c) 


Nhf 


a half of 5*Ni amount: Mm = O.OISM© 


Nno 


no radioactive '*Ni 


Nvm 


'*Ni mixed up to a velocity of 830 km s~' 


Nvh 


"'''Ni mixed up to a velocity of 1095 km s"* 


Lop 


no contribution of the expansion opacity 


Ldk 


without the limb-darkening effect 



is in equilibrium with matter, the radiative diffusion reduces 
both the radiation and internal gas energies. By day 116 the 
total optical depth of the ejecta falls to a value of ~ 4, and 
a low opacity of matter significantiy weakens the interaction 
between radiation and matter. As a consequence, for the next 
9 days the photon diffusion mainly reduces the stored radiation 
energy and eventually exhausts it, and the luminosity decreases 
abruptly and becomes close to the gamma-ray luminosity due 
to the decay of the radioactive ^*Co, which is calculated by tak- 
ing account of the retardation effect in the envelope transparent 
for gamma-rays (Fig. 13a). 

5.7. Plateau tail phase 

By the end of the previous phase, the bolometric luminosity 
does not fall directly to the gamma-ray luminosity (Fig. 13a), 
and the total radiation energy is not exhausted completely. 
Nearly at the same time the photosphere disappears, and the 
SN envelope becomes optically thin for radiation. For exam- 
ple, the total optical depth of the ejecta is 2.96 at the moment 
ts, while the optical depth above the ^^Ni bubble shell is only 
0.78. Because during the phase of the radiation energy exhaus- 
tion the photon gas in the outer layers cools more intensively 
than in the inner part of the envelope (t = 126 and 130 days 
in Fig. 13b), the radiation flow is generated in the warm imier 
layers (t = 130 and 134 days in Fig. 13b), then it propagates 
throughout the transparent layers, and results in some luminos- 
ity excess in the light curve just after the steep decUne in the 
luminosity compared to the radioactive tail (Fig. 13a). We call 
this behavior of the tight curve after the main plateau at the 
beginning of the radioactive tail "a plateau tail". 

The radiation flow generated near the center of the envelope 
exists from f ~ 130 days to t ~ 150 days and then disappears 
(Fig. 13b). The characteristic duration of the plateau tail in the 
optimal model is 26 days (Fig. 13a). Note that the energy ex- 
cess radiated during the plateau tail over the luminosity of ra- 
dioactive decays is 2.5 x 10"*^ erg and 1 1% of the total radiation 
energy within the ejecta at the moment ts . 



5.8. Radioactive tail of the light curve 

After the phase of the plateau tail, the bolometric luminosity 
of the optimal model decreases to the gamma-ray luminosity 
due to the decay of the radioactive ^^Co (Fig. 13a) and is solely 
powered by the radioactive energy source. This is how the stage 
of the radioactive tail starts, and the calculated light curve is in 
good agreement with the observations of SN 1999em after day 
150 (Fig. 3). Note that the envelope remains optically thick for 
the gamma rays during a few hundred days (the total optical 
depth is 30.6 at t = 150 days), and they deposit the bulk of 
their energy locally, at least at the beginning of the radioactive 
tail. 

6. General properties of the optimal model 

It is well known that the photometric characteristics of the 
SN IIP outburst are determined mainly by the basic parameters: 
the pre-SN radius Rq, the ejecta mass M,>„, , and the explosion 
energy E. Moreover, to get better agreement with the observa- 
tions of SN 1999em, we have to vary the density distribution 
in the pre-SN model, its chemical composition, the ^^Ni mass, 
and its mixing. A dependence of the bolometric light curve on 
the expansion opacity and the effect of limb darkening is of par- 
ticular interest. These general properties of the optimal model 
are studied by means of the models listed in Table 2. 

6.1. Presupernova structure 

At first, we study the dependence of the bolometric luminosity 
on the inner layers' density in the pre-SN model and consider 
the extreme case of an auxiliary model similar to the optimal 
model but without a dense central core (Fig. 14a). In such a 
model the inner layers do less work moving outwards in the 
gravitational field than do those in the optimal model with a 
dense central core and hence acquire higher velocities than in 
the optimal model. In contrast, the outer layers of the optimal 
model expand faster than those of the auxiliary model. As a 
consequence, the photospheric radius of model Dll is greater 



10 



V. p. Utrobin: Type IIP supernova 1999em 




Fig. 14. Dependence on the pre-SN structure. Top: (a) density distribu- 
tion with respect to interior mass for the pre-SN model Dl 1 (solid line) 
and the similar model Dec but without a dense central core (dashed 
line), and (b) the corresponding bolometric light curves. Bottom: (c) 
density distribution with respect to radius for the pre-SN model Dll 
(solid line) and the similar models Ldn and Mdn, but with less (dotted 
line) and more (dashed line) dense outer layers, respectively, and (d) 
the corresponding bolometric light curves. 



than in the auxiliary model up to about day 35 and is smaller 
later on. At the comparable effective temperatures of the mod- 
els, this results in completely different behavior by the bolo- 
metric light curves (Fig. 14b). This implies that there is a dense 
central core inside the real pre-SN of SN 1999em. 

In terms of the CRW properties the growth in the bolomet- 
ric luminosity during the CRW phase in the auxiliary model 
(Fig. 14b) is controlled by effective index n at the photosphere 
position that is greater by ~ 3.5-0.5 than in the optimal model. 
A greater effective index corresponds to a greater increase in 
the total flux of mass through the surface of the CRW front with 
time and, as a consequence, to faster growth in the bolometric 
luminosity with time (Grassberg & Nadyozhin 1976). 

The influence of the density distribution in the outer lay- 
ers of the pre-SN on the bolometric luminosity is illustrated 
by model Dll and the same models but with both less dense 
and denser outer layers (Figs. 14c and 14d). A roughly simi- 
lar behavior of the effective index n at the photosphere position 
in these models reflects the approximately equal slopes of the 
corresponding bolometric light curves during the CRW phase. 
A transition from the model with the less dense outer layers 
to that of the denser layers is characterized by a reduction in 
the cooling by the adiabatic expansion and, consequently, by 
an increase in the bolometric luminosity in the whole SN out- 
burst except for phases of the shock breakout and the radioac- 
tive tail. As a result, for these models the energies radiated dur- 
ing the first 180 days are 1.251 x lO'^' erg, 1.355 x erg, and 



Fig. 15. Dependence on the chemical composition in the outer layers 
of the ejecta beyond the helium core. Top: (a) the chemical composi- 
tion of model Hoi that differs from model Dl 1 in the mass fraction of 
hydrogen (solid line) and helium (long dashed line) (Table 2), the dis- 
tribution of heavy elements (short dashed line) and radioactive '''Ni 
(dotted line) being the same, and (b) the corresponding bolometric 
light curve (dotted line) and that of model Dll (solid line). Bottom: 

(c) the chemical composition of model Zol that differs from model 
Dl 1 by the mass fraction of helium and heavy elements (Table 2), and 

(d) the corresponding bolometric light curve (dotted line) and that of 
model Dll (solid line). 



1.430 X 10"*^ erg, respectively. Thus, the bolometric Ught curve 
is fairly sensitive to the initial structure of the outer layers. 

6.2. Chemical composition 

In addition to the pre-SN structure, the chemical composition 
of the ejecta also affects the SN luminosity. Let us consider a 
dependence on the chemical composition in the outer layers of 
the ejecta beyond the heUum core. An enhancement of helium 
abundance at the expense of hydrogen abundance in model Hoi 
(Fig. 15a and Table 2) results in the higher luminosity during 
the adiabatic cooling phase than that of model Dll, nearly the 
same luminosity in the CRW phase, a small bump at the end 
of the plateau and then the shortened duration of the plateau 
(Fig. 15b). This behavior of the luminosity is due to the lower 
opacity in the outer layers of model Hoi. On the other hand, 
lower metallicity for model Zol (Fig. 15c and Table 2) than 
for model Dll (Table 1) increases the bolometric luminosity 
slightly during the CRW phase and then reduces it somewhat 
at the phase of radiative diffusion cooling (Fig. 15d). This is ex- 
plained by a smaller contribution of heavy elements to opacity 
in the outer layers of the ejecta in model Zol. Thus, the bolo- 
metric light curve depends weakly on the chemical composi- 
tion in the outer layers of the ejecta beyond the helium core. 
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Fig. 16. Dependence on the chemical composition in the inner layers 
of the ejecta. Top: (a) the chemical composition of model Sci, which 
differs from model Dll by the sharp metals/He/H composition inter- 
face, and (b) the corresponding bolometric light curve (dotted line) 
and that of model Dll (solid line). Bottom: (c) the chemical com- 
position of model Hec, which differs from model Dll by the larger 
helium core, and (d) the corresponding bolometric light curve (dotted 
line) and that of model Dll (solid line). 



Fig. 17. Influence of the '*Ni mass and its mixing. Panel (a): bolo- 
metric light curves for model Dll (solid line) and the similar models 
Nhf and Nno but calculated with a half of the ^''Ni amount (dotted 
line) and without '*Ni (dashed line), respectively. Panel (b): bolomet- 
ric light curves at the end of the plateau for model Dll (solid line), 
model Nvm (dashed line), and model Nvh (dotted line) in which ^*Ni 
is mixed up to a velocity of 660, 830, and 1095 kms"', respectively. 



ing the evolutionary state of a massive star. In this context it is 
intriguing that the increase in the helium core mass from 5.6 
Mq in model Dll to 8.1 M© in model Hec (Fig. 16c), under 
a deep hydrogen mixing downward to x 700 kms"', leads to 
some decrease in the opacity of matter in the mixing region 
at the heUum/hydrogen composition interface and, as a con- 
sequence, to the increase in the bolometric luminosity at the 
plateau (Fig. 16d). The resulting bolometric light curve is still 
consistent with the observations of SN 1999em. Thus, the mass 
of the helium core between 5.6 Mq and 8.1 Mq is acceptable 
for the optimal model. 



while hydrogen is abundant there and controls opacity of mat- 
ter. 

A good fit of the calculated bolometric light curve for the 

optimal model to what is observed (Fig. 3) is obtained with mu- 
tual mixing of the hydrogen-rich and helium-rich matter in the 
inner layers of the ejecta (Fig. 2). It is appropriate to ask here 
in what way the sharp boundary between the hydrogen-rich and 
helium-rich layers at the edge of the helium core — a character- 
istic of the evolutionary models of pre-SNe — changes the Ught 
curve. In the auxiliary model Sci without the mutual mixing at 
the edge of the helium core (Table 2) this boundary is located at 
the interior mass of 5.4 Mq (Fig. 16a). It is known that a transi- 
tion from the hydrogen-rich to helium-rich layers, i.e. from the 
low to high ionization potential matter, causes an increase in 
the effective temperature and, as a consequence, a growth in the 
bolometric luminosity (Grassberg & Nadyozhin 1976; Utrobin 
1989). In fact, when the CRW front reaches the sharp bound- 
ary between the hydrogen-rich and helium-rich layers, crosses 
it, and then enters the helium core, a distinct bump in the light 
curve at the end of the plateau appears (Fig. 16b). Evidently, 
this local maximum of the luminosity is inconsistent with the 
observations of SN 1999em and other SNe IIP. 

It is well known that at the time of an SN event the lu- 
minosity of a massive progenitor is determined mainly by the 
mass of the helium core, since a contribution of the hydrogen 
burning shell to the total energy generation is negligible. Hence 
the mass of the heUum core is a critical quantity in evaluat- 



6.3. Mass of^^Ni and its mixing 

After the CRW stage, when the radiative difl'usion takes place, 
the radioactive decay of the ^''Ni and ^^Co nuclides begins 
to dominate in powering the luminosity. This fact is clearly 
demonstrated by model Nno without ^^Ni in the envelope 
(Fig. 17a and Table 2). A lack of ^^Ni shortens the duration 
of the bolometric light curve to about 100 days. It is evident 
that at later times the bolometric light curve depends on both 
the total mass of ^^Ni and its distribution over the ejecta. Of 
course, the total mass of ^^Ni is measured by the radioactive 
tail of the observed Ught curve and is 0.036 Mq in the optimal 
model for the Cepheid distance of 1 1 .7 Mpc. A reduction in the 
total mass of ^''Ni from 0.036 Mq (model Dll in Table 1) to 
0.018 Mq (model Nhf in Table 2) results mainly in a shorter 
duration of the light curve (Fig. 17a). Thus, the duration of the 
bolometric light curve has a minimum value of about 100 days 
and is a function of the ^*Ni amount. 

To investigate the influence of the ^^Ni mixing on the light 
curve, hydrodynamic models Nvm and Nvh were calculated 
with the ^*Ni distribution difl'erent from that of model Dll, but 
containing the same total mass of ^''Ni (Table 2). In these mod- 
els the radioactive nickel is spread over a larger mass range 
of the ejected envelope or, equivalently, over a wider velocity 
range than in model Dll. Hence, more favorable conditions for 
the radiative diffusion and an additional heating due to the ra- 
dioactive decays should appear and therefore should increase 
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Table 3. Basic parameters at the phase of the plateau tail. 



Model 














AE^//E'f 




(days) (days) (days) 






(lO"* erj 


I) 


Dll 


124 


26 


8.76 


2.97 


2.96 


2.5 


0.110 


Lop 


110 


25 


7.83 


3.19 


1.71 


3.8 


0.136 


Nhf 


115 


22 


7.94 


2.77 


2.99 


1.5 


0.100 


Mps 


127 


22 


7.89 


2.79 


2.64 


2.3 


0.089 


Bps 


119 


28 


9.25 


3.03 


3.01 


3.1 


0.107 



Fig. 18. Dependence of the plateau tail of the bolometric light curve 
on (a) the '*Ni amount (model Nhf in Table 2, dotted line) and the 
contribution of the expansion opacity (model Lop in Table 2, dashed 
line), and on (b) the ejecta mass (model Mps in Table 4, dotted line) 
and the explosion energy (model Hps in Table 4, dashed line). Thick 
solid line is the bolometric light curve of model Dll and thin solid 
line is the luminosity of radioactive decays. 

the luminosity at the phase of radiative diffusion cooling that 
is clearly seen in Fig. 17b. In addition, this increase in the lu- 
minosity causes its subsequently earlier fall to the radioactive 
tail. As a result, with increasing the ^*Ni mixing in velocity 
space the linear luminosity shoulder at the phase of radiative 
diffusion cooling transforms into the convex light curve at the 
end of the plateau. Thus, a linear shoulder of the luminosity in 
logarithmic scale at the phase of radiative diffusion cooling im- 
plies a weak ^*Ni mixing. Because a discrepancy between the 
light curve of model Dll (which is in agreement with the ob- 
servations of SN 1999em) and those of models Nvm and Nvh 
is rather significant, one can conclude that the bulk of the ra- 
dioactive ^^Ni should be confined to the innermost layers of 
the ejected envelope expanding with velocities less than 660 
kms"^ 

6.4. Plateau tail 

The general properties of the plateau tail on the radioactive tail 
are illustrated by the light curves shown in Fig. 18 and the rel- 
evant parameters of the models listed in Table 3 where is the 
time of the onset of the plateau tail as defined in Fig. 6, Af^'' 
the characteristic duration of the plateau tail, AIr - Rext/c is 
the characteristic time of radiation flow for the external radius 
of the expelled envelope Rgxi in the middle of the plateau tail, 
Tfof the total optical depth of the ejecta at the moment 1$, AEj'' 
the energy excess radiated during the plateau tail over the lu- 
minosity of radioactive decays, and AE^''/E'°' is a ratio of this 
energy to the total radiation energy within the ejecta at the mo- 
ment ts. 

It is clear that there is no correlation between the character- 
istic duration of the plateau tail and the total optical depth of 
the ejecta (Table 3), and the duration depends very weakly on 
the opacity of matter (models Dll and Lop in Table 3). These 
regularities reflect the fact that the radiation responsible for the 
plateau tail flows throughout an optically thin medium. On the 
other hand, the characteristic duration is proportional to the 
characteristic time of radiation flow independent of the param- 
eters of the listed models (Table 3), and their ratio is 2.95 on 
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Fig. 19. The bolometric Ught curves for model Dll (solid line) and 
the similar models Lop and Ldk, but calculated without the expansion 
opacity (dashed line), and for the isotropic emergent radiation, without 
the Umb darkening (dotted line). 

average. It is worth noting that the energy excess of the plateau 
tail is roughly 10% of the total radiation energy at the moment 
ts for all but model Lop (Table 3). 

The plateau tail of the bolometric light curve demonstrates 
the following important properties. A decrease in the ^^Ni mass 
(models Dll and Nhf) shortens the characteristic duration of 
the plateau tail and reduces its energy excess (Fig. 18a and 
Table 3). This influence of the ^^Ni amount is related to its large 
role in powering the light curve at the end of the main plateau 
after day 100 (Fig. 17a). The higher the average expansion ve- 
locity of the ejecta (models Mps, Dll, and Bps), the longer 
the duration of the plateau tail and the higher the energy ex- 
cess (Fig. 18b and Table 3). An increase in the duration of the 
plateau tail with the average expansion velocity results from 
its proportionality to the characteristic time of radiation flow. 
The higher average expansion velocity causes the onset of the 
plateau tail at an earlier time, when the total radiation energy is 
higher and, as a consequence, the energy excess, radiated over 
the luminosity of radioactive decays during the plateau tail, be- 
comes greater. 

6.5. Expansion opacity 

A contribution of numerous metal lines to opacity plays a fun- 
damental role in reproducing the observed light curve of SN 
1999em, a normal SN IIP, as shown by the optimal model Dll 
and model Lop calculated by neglecting the expansion opac- 
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Fig. 20. Evolution of the velocity distribution in mass for models Dl 1 
{solid line) and Lop (dotted line): at f = 0.6234 days when the shock 
approaches the pre-SN surface, at f = 0.8937 days when the shock 
emerges on the surface, at t = I days, at f = 1.5 days, and sub- 
sequently, at later times t > 10 days when the envelope is already 
expanding homologously. The mass is measured from the envelope 
surface. 



ity (Fig. 19). The decrease in opacity due to disregarding the 
line contribution in model Lop speeds up the radiation diffusion 
and causes the bolometric luminosity to increase significantly 
compared to that of model Dl 1 during nearly the whole out- 
burst and, as a consequence, the characteristic duration of the 
plateau to shorten by about 10 days. It is equivalent in action to 
an increase in the explosion energy that will be demonstrated 
in Sect. 7 (Fig. 22c). This fact leads to the firm conclusion that 
neglecting the expansion opacity underestimates the explosion 
energy and introduces an error of nearly 20% into its value. 

It is worth noting that the slope of the light curve shoulder 
at the phase of radiative diffusion cooling in model Dl 1 is flat- 
ter than in model Lop without the expansion opacity (Fig. 19). 
This behavior of the light curve is consistent with the depen- 
dence of the characteristic diffusion time on the optical depth 
of the diffusion region (1) which, in turn, is proportional to the 
opacity of matter 

The resonance scattering of radiation in numerous metal 
Unes is also important in forming the gas flow in the outer lay- 
ers of the expanding SN envelope (Fig. 20). After the shock 
emerges on the pre-SN surface at f = 0.8937 days, an expan- 
sion of the envelope begins, and it becomes homologous during 
the next 10 days. An additional acceleration induced by the ex- 
pansion opacity occurs between f » 1 days and f 1.5 days, 
lasts up to ? « 7.5 days, and takes place in the outer layers 
with a mass of about 10""* Mq. A comparison of models Dll 
and Lop shows that these layers gain a velocity excess as large 
as 600 kms"\ except for the outermost layers with a mass of 




time 



Fig. 21. Schematic bolometric light curve of SN IIP. Lettered dots 
mark the onset of the luminosity rising (o), the middle of the plateau 
(p), and the end of the plateau (e). 

about 10"** Mq. Note that this additional acceleration is respon- 
sible for a feature in the gas flow as seen in Fig. 9 where the 
velocities of different mass shells in the ejected envelope as a 
function of time are overplotted on the photospheric velocity. 

6.6. Effect of limb darkening 

During the first 116 days, when a well-defined photosphere 
exists, particularly at the CRW phase, the emitted radiation 
is nearly isotropic, and the bolometric light curves of models 
Dll and Ldk, calculated without the effect of the limb dark- 
ening, almost coincide (Fig. 19). As the envelope expands and 
becomes optically thin, the continuum formation region grad- 
ually becomes more extended, and the degree of anisotropy of 
the emergent radiation increases. The increase in the radiation 
anisotropy with time clearly results in the difference between 
the hght curves of models Dll and Ldk during the phase of 
the radiation energy exhaustion: the bolometric luminosity cal- 
culated by taking the limb-darkening law into account is lower 
than for isotropic radiation (Fig. 19). 

Utrobin (2004) shows that the difference between the bolo- 
metric luminosity calculated by taking into account the emer- 
gent anisotropic radiation and the luminosity for the isotropic 
radiation depends on both the limb-darkening law and the retar- 
dation effect. Interestingly, this difference increases with grow- 
ing the degree of the anisotropy of the emergent radiation, 
while its sign is determined solely by the time derivative of 
the bolometric luminosity. The latter property explains the be- 
havior of the bolometric luminosity for models Dll and Ldk 
discussed above. It seems that in SN 1999em, a normal SN IIP, 
the effect of limb darkening is more of purely academic interest 
than of practical value. 

7. Physical and observed parameters 

Generally, our main goal is to construct the adequate hydro- 
dynamic model by fitting the photometric and spectroscopic 
observations of the object under study. Sometimes it is tempt- 
ing to evaluate the relevant physical parameters for other SN 
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with close observed properties. In this case, it is reasonable 
to use the results already obtained for the similar object. For 
this purpose, we have to find the relationships between the ba- 
sic physical and observed parameters. Evidently, the important 
physical parameters are the initial radius R(}, the ejecta mass 
Menv, and the explosion energy E. The large impact of ^^Ni on 
the Ught curve near the end of the plateau, which is shown and 
discussed in Sect. 6.3, makes us take the total ^''Ni mass into 
account. Because of the temporal character of the energy de- 
position from radioactive decays, it is impossible to treat this 
energy input as an addition to the explosion energy, so we have 
to introduce the total ^^Ni mass M^i as the fourth basic physical 
parameter. 

To describe the observed properties of SN IIP, Litvinova & 
Nadyozhin (1983, 1985) used the three parameters: the plateau 
duration, the absolute V magnitude, and the photospheric ve- 
locity in the middle of the plateau. However, in the course of 
time it became evident that this simple description needed to 
be corrected for the influence of ^^Ni on the light curve via 
an additional parameter (Nadyozhin 2003). Unfortunately, the 
bounds of the plateau are unclear because they are no specific 
points on the Ught curve. Moreover, the interval of a plateau 
duration misses the bolometric luminosity jump at the shock 
breakout and the subsequent adiabatic cooling phase together 
with the relevant dependence on the initial radius, the ejecta 
mass, and the explosion energy (Fig. 22). To overcome this 
problem, instead of the plateau duration we define a charac- 
teristic duration At of the light curve as a whole, measuring it 
from the bolometric luminosity jump at the shock breakout to 
the end of the plateau for an e-fold reduction of the luminos- 
ity in the middle of plateau (Fig. 21). Of course, in this way 
the characteristic duration is easy to determine for a theoreti- 
cal hght curve, as shown in a schematic bolometric light curve 
of SN IIP, but very difficult for an observed one. Fortunately, 
this situation is not hopeless. First, it is possible to apply the 
theoretical light curve as a template for measuring the specific 
moments to and te (Fig. 21) and, consequently, for calculating 
the characteristic duration. Then, as the observational data be- 
come more and more extensive, the relevant templates may be 
constructed from the observed light curves of the well-studied 
SNe IIP. The second modification we made is to use the bolo- 
metric luminosity L^^^ in the middle of the plateau because it 
is more adequate than the absolute V magnitude for estimating 
the explosion energy. The photospheric velocity v|J^^ in the mid- 
dle of the plateau is the third observed parameter in accordance 
with Litvinova & Nadyozhin (1983, 1985). 

To establish the relationships among the basic physical and 
observed parameters accepted above, we carried out a param- 
eter study of hydrodynamic models varying the physical pa- 
rameters and measuring the observed ones. Evidently, a four- 
parameter approximation is a crude description of the hydro- 
dynamic model properties. The sensitivity of the bolometric 
light curve to the pre-SN structure and its chemical compo- 
sition makes this approximation valid only for the vicinity of 
the optimal model Dl 1, not for a whole region of the basic pa- 
rameters. We restricted ourselves to a local parameter study in 
the vicinity of the optimal model and explored a 15% limited 
range of the basic physical parameters (Table 4). The result- 



Table 4. Physical parameters of auxiliary hydrodynamic models. 



Model 


Ro 




E 




X 


Z 




(Ro) 


(Mq) (10" eq 


;) (lO-^Mo) 






Dll 


500 


19 


1.30 


3.60 


0.735 


0.017 


Rms 


425 


19 


1.30 


3.60 


0.735 


0.017 


Rps 


575 


19 


1.30 


3.60 


0.735 


0.017 


Mms 


500 


16 


1.30 


3.60 


0.735 


0.017 


Mps 


500 


22 


1.30 


3.60 


0.735 


0.017 


Ems 


500 


19 


1.10 


3.60 


0.735 


0.017 


Eps 


500 


19 


1.50 


3.60 


0.735 


0.017 


Nms 


500 


19 


1.30 


3.06 


0.735 


0.017 


Nps 


500 


19 


1.30 


4.14 


0.735 


0.017 
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Fig. 22. Dependence of the bolometric light curve of the optimal 
model Dll (solid line) on the basic parameters: (a) the initial radius, 
model Rms {dotted line) and model Rps {dashed line); (b) the ejecta 
mass, model Mms {dotted line) and model Mps {dashed line); (c) the 
explosion energy, model Ems {dotted line) and model Eps {dashed 
line); (d) the total '*Ni mass, model Nms {dotted line) and model Nps 
{dashed line). 



ing bolometric light curves of auxiliary hydrodynamic models 
are plotted in Fig. 22 and compared to the optimal model. The 
basic observed parameters of these models are evaluated ac- 
cording to the definition, given in Fig. 21, and listed in Table 5. 

The influence of the basic physical parameters is illustrated 
by the bolometric light curves of the auxiliary models (Fig. 22). 
Increasing the initial radius causes an increase in the character- 
istic time of the envelope expansion and, as a consequence, an 
increase in the width of the narrow intense peak on the bolo- 
metric light curve produced by the heating of the outer layers 
at the shock breakout phase. In addition, an increase in the ini- 
tial radius leads to a reduction in the cooling by the adiabatic 
expansion and, accordingly, to an increase in the bolometric lu- 
minosity during the whole SN outburst except for the phase of 
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Table 5. Observed parameters of auxiliary hydrodynamic models. 



log£ = -1.4761ogL^^, + 4.8991ogv^^+3.0401ogA? 
+0.3121ogMNi + 90.765, 



(4) 



Model 


A f 




ph 


Ikin 


Irad 




(days) 


(ergs"') 


(kms"') 






Dll 


117.57 


42.0350 


2374.4 


0.8541 


0.8258 


Rms 


114.88 


41.9927 


2335.3 


0.8088 


0.8429 


Rps 


118.25 


42.0521 


2388.3 


0.8354 


0.7826 


Mms 


111.70 


42.0542 


2481.5 




0.7881 


Mps 


121.64 


42.0170 


2291.9 


0.9141 


0.8252 


Ems 


122.50 


41.9629 


2134.7 


0.8134 


0.8040 


Eps 


111.93 


42.0979 


2640.5 


0.8890 


0.8424 


Nms 


113.98 


42.0372 


2452.5 


0.8943 


0.8303 


Nps 


119.39 


42.0264 


2320.3 


0.7968 


0.8102 



the radioactive tail (Fig. 22a). It is interesting that this behav- 
ior by the bolometric light curve with initial radius is similar 
to its dependence on the density structure in the outer layers 
(Fig. 14c), which is discussed in Sect. 6.1. The latter depen- 
dence may be considered as a dependence on some effective 
radius of the pre-SN. In other words, an increase in the density 
in the outer layers of the pre-SN mimics an increase in the ini- 
tial radius. Both decreasing the ejecta mass and increasing the 
explosion energy enlarge the average expansion velocity of the 
envelope. The increase in the average velocity of the envelope 
makes the shock wave propagating through the pre-SN matter 
stronger and, as a result, it heats the matter to a higher temper- 
ature, which increases the luminosity in the narrow peak of the 
light curve, at the subsequent phase of adiabatic cooling and at 
the beginning of the CRW phase (Figs. 22b and 22c). By the 
end of the CRW phase, decreasing the ejecta mass reduces the 
bolometric luminosity and shortens the characteristic duration 
of the hght curve (Fig. 22b). In contrast, increasing the explo- 
sion energy results in the luminosity growing and, accordingly, 
in shortening the characteristic duration as well (Fig. 22c). Up 
to the phase of the radiative diffusion cooling, the radioactive 
^^Ni does not affect the light curve; but during this phase and 
later on, the energy deposition from radioactive decays pow- 
ers the luminosity (Fig. 22d). The higher the total ^^Ni mass, 
the greater the energy deposition and the later the luminosity 
descending to the radioactive tail. Note that a slope of the lu- 
minosity shoulder at the phase of the radiative diffusion coohng 
clearly depends on the initial radius (Fig. 22a) and on the ex- 
plosion energy (Fig. 22c) in a qualitative agreement with the 
simple description of photon diffusion (1). 

To evaluate the three basic physical parameters (the initial 
radius, the ejecta mass, and the explosion energy) from the four 
values measured from observations (the characteristic duration 
of the light curve, the bolometric luminosity and the photo- 
spheric velocity in the middle of the plateau, and the total ^^Ni 
mass), we use the following approximate relations: 

Iog/?o = +3.4811ogL{;^^,- 5.937 log v^^- 1.999 log A? 

-0.499 logMNi - 120.174 , (2) 

log Men. = -2.942 log Ll^i + 7.606 log v^^ +7.807 log ht 

-0.042 log Mni + 83.045 , (3) 



where Rq, Menv, and Mni are in solar units, E in units of erg, L^^^ 
in erg s"', v|J^^ inkms"', and Af in days. The above approximate 
formulae (2-4) are found by means of minimizing the errors 
between the values calculated with these relations and the cor- 
responding physical parameters from Table 4. For this purpose 
we constructed a variational functional in a quadratic form de- 
pending on the functional values, i.e. coefficients in the approx- 
imate formulae (2-4), and then used the direct search method 
of Powell ( 1964). Note that the obtained approximate formulae 
give the estimates of the physical parameters for the auxiliary 
models, the accurate values of which are listed in Table 4, with 
an error as great as » 2.5%. 

In addition to the observed parameters of the auxiliary mod- 
els, the important energetic ratios qkm - 0.5Me„i.v^^^ /£■,(:,„ and 
^rad = L^gi^t/Erad, where Ellin is the kinetic energy of the 
ejecta and Erad is the total radiation energy emitted during the 
first 180 days, are given in Table 5. The mean values of the 
ratios qidn and qrad are 0.842 and 0.816 with the standard de- 
viations of 4.4% and 2.3%, respectively. A quite low value of 
the standard deviation in the energetic ratios shows that the ki- 
netic energy of the envelope and the total radiation energy are 
approximated more or less accurately by the observed parame- 
and Af. 

A functional dependence of the observed parameters on the 
physical ones is, in turn, derived directly from the approximate 
formulae (2-4) and can be written as 



ters: if „ v^, 

boP ph 



logLl^i = -1-0.488 log/?o - 0.267 log Me„v + 1.006 log £ 
-0.0821ogMNi- 10.469, 

logv^^ = +0.083 log^o- 0.247 log M,, 
-0.1841ogMNi- 31.951 , 



•0.688 log £ 



logAf = +0.103 log/?o + 0.268 log M^, 
+0.1541ogMNi+ 16.548. 



0.291 log £ 



(5) 



(6) 



(7) 



The estimates of the observed parameters for the auxiliary 
models obtained with the above relations (5-7) differ from the 
accurate values given in Table 5 by an error less than « 2.2%. 

Finally, it should be emphasized that all of the above for- 
mulae approximate the physical parameters and the observable 
properties of hydrodynamic models but only in the vicinity of 
the optimal model Dll. To cover a wider range of the parame- 
ters, we have to construct a global approximation, which should 
be based, in our opinion, on a set of hydrodynamic models of 
many real SNe IIP, but not on that of arbitrary hydrodynamic 
models. 

8. Comparison with SN 1987A 

Now we have got a confidence in reproducing the observed 
bolometric light curve and the Ha line of SN 1999em, a nor- 
mal SN IIP, we can compare this object with the well-studied 
SN 1987A, a pecuUar SN IIP. The important parameters of the 
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Table 6. Hydrodynamic models for SN 1999em and SN 1987 A. Table 7. Shock breakout and UV flash in SN 1999em and SN 1987 A. 



SN Ro M,„„ E Mn, Z vj;f^ v™« 

(Rq) (Mq) (10^' erg) (10-2 M©) (kms-') (kms-') 



99eni 500 19 1.3 3.60 0.017 660 700 
87A 35 18 1.5 7.65 0.006 3000 600 



SN t,h T^;; logL^," Ml N"" f E^^ 

(days) (105 K) (10^ K) (erg s"') (days) (lO^*) (days) (lO^* erg) 

99em 0.865 3.84 1.76 44.81 0.0200 276.8 1.234 159 
87A 0.049 13.2 5.68 44.63 0.0007 7.615 0.266 6.24 



optimal model Dll for SN 1999em and model M18 for SN 

1987A (Utrobin 2005) are given in Table 6. It is remarkable 
that the ejecta mass Menv, the explosion energy E, and the min- 
imum velocity of the hydrogen-rich envelope v™" are compa- 
rable, while the pre-SN radius Rq, the mass fraction of heavy 
elements in the outer layers of the pre-SN Z, the total mass of 
^^Ni Mm, and its maximum velocity v™* are quite different. 

The mass of a star is a fundamental parameter that deter- 
mines the properties of the star and the course of its evolution. 
The pre-SN masses of SN 1999emandSN 1987A, estimated by 
the mass of the neutron star and the ejecta mass, are nearly 20.6 
Mq and 19.6 Mq, respectively. The masses of helium cores in 
the pre-SNe of SN 1999em and SN 1987A are 5.6-8.1 Mq and 
6.0 Mq (Woosley 1988), respectively. These masses are close 
enough to suppose that, in the final stages of stellar evolution, 
nearly the same iron cores form within the pre-SNe. This fact 
and roughly the same explosion energies of SN 1999em and 
SN 1987A imply a unique explosion mechanism for these core 
collapse SNe. 

It is evident that the amount of radioactive ^*'Ni and its dis- 
tribution in the SN envelope are a clear trace left by the explo- 
sion mechanism. Both SN 1999em and SN 1987A differ in the 
total mass of radioactive ^^Ni and in its distribution throughout 
the envelope (Table 6). Note that for these SNe the total ^*'Ni 
mass is correlated to the explosion energy. This dependence is 
consistent with the empirical correlation between the total ^^Ni 
mass and the explosion energy for SNe IIP found by Nadyozhin 
(2003) and Hamuy (2003). The optimal model for SN 1999em 
is characterized by a weak ^^Ni mixing in velocity space up to 
^ 660 km s ' and a deep hydrogen mixing downward to ^ 700 
kms"' (Fig. 2 and Table 6). In turn, the hydrodynamic model 
for SN 1987 A is characterized by a moderate ^^Ni mixing up 
to ~ 3000 kms ' and a deep hydrogen mixing downward to 
~ 600 kms-i (Uti-obin 2005). 

The structure of the pre-SN and the chemical composition 
of its outer layers, which are unaffected by explosive nucle- 
osynthesis and which are the end result of the entire evolution 
of the star, determine the pattern of SN outburst in many re- 
spects. The pre-SN model for SN 1999em has the structure of a 
red supergiant with a radius of 500 Rq, and the chemical com- 
position of its outer layers is the standard solar composition. 
It is quite definite that the progenitor of SN 1987A was the 
star Sanduleak -69°202, a B3 la blue supergiant. The pre-SN 
model for SN 1987 A has a radius of 35 Rq, and the chemical 
composition of its outer layers is typical of the LMC chemical 
composition, X = 0.743, Y = 0.251, and Z = 0.006 (Dufour 
1984). It is worth noting that a deficit of heavy elements in the 
LMC matter compared to the standard solar composition favors 
the formation of blue supergiants (Amett 1987; Hillebrandt et 



al. 1987). In the case of SN 1987A, a relative compactness of 
the pre-SN is a major factor in understanding the pecuUar prop- 
erties of this phenomenon (Grassberg et al. 1987). 

The basic characteristics of the shock breakout and the UV 
flash in SN 1999em and SN 1987A are summarized in Table 7 
where tsh is the time when the shock wave reaches the stellar 
surface, the maximum of the color temperature, the 
maximum of the effective temperature, L™"^^ the maximum of 
the bolometric luminosity, AIl the width of the luminosity peak 
at a half level of its maximum, N"" the total number of ionizing 
photons above 13.598 eV for the whole outburst, t'^ the time 
when the number of ionizing photons measures up 90% of the 
total number, and the energy radiated by the SN during 
this time. 

Under the comparable ejecta masses and the explosion en- 
ergies, the smaller the pre-SN radius, the higher the velocity of 
matter in the outermost layers with a sharp decline of density 
where the matter accelerates due to the effect of hydrodynamic 
cumulation. At the same time the characteristic time of the en- 
velope expansion becomes shorter and the adiabatic losses of 
energy greater. Clearly, in SN 1999em the shock wave reaches 
the stellar surface much later than in SN 1987A, and the color 
and effective temperatures jump to lower values, while the 
bolometric luminosity rises to a higher value and its peak is 
much wider (Table 7). Note that in SN 1999em the bolometric 
luminosity peak coincides with the maximum of the effective 
temperature and occurs after the shock wave reaches the stellar 
surface in contrast to SN 1987A in which the maximum of the 
effective temperature precedes the peak, and the shock wave 
emerges between them on the pre-SN surface. The UV flash in 
SN 1999em produces a greater total number of ionizing pho- 
tons by a factor of 36.3, and it lasts longer than in SN 1987A 
by a factor of 4.6, the radiated energy being greater by a factor 
of25.5(Table7). 

The smaller the pre-SN radius, the greater the adiabatic 
losses of energy and the shorter the CRW phase — a charac- 
teristic feature of SNe IIP. A difference between the explosions 
of the red and blue supergiants is radical in the total radiation 
and gas energies as a function of time (Fig. 23) and, especially, 
in the Ught curves (Fig. 24). To clarify the physical processes 
in SN 1999em and SN 1987 A, we compare the optimal models 
of these SNe with the special hydrodynamic models containing 
no ^^Ni in the envelope. Note that during the entire outburst the 
total radiation energy is much greater than the total gas energy 
of the envelope. A lack of ^^Ni makes the energy losses due 
to the adiabatic expansion a dominant process controlling the 
energy balance of an SN. In this case, the blue supergiant struc- 
ture of the pre-SN for SN 1987A results in such huge energy 
losses that its total radiation energy is much less than that of SN 
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Fig. 23. Total radiation (thick solid line) and gas (thin solid line) en- 
ergies as a function of time for SN 1999em and SN 1987A. The de- 
pendence on the energy deposition from radioactive decays is shown 
by the special hydrodynamic models without ^*Ni in the envelope 
(dashed line). 



Fig. 24. Comparison of the calculated bolometric light curve (thick 
solid line) with the gamma-ray luminosity (thin solid line) for SN 
1999em and SN 1987 A. The special hydrodynamic models contain- 
ing no '*Ni in the envelope show the influence of the energy deposition 
from radioactive decays (dashed line). 



1999em the whole time and is completely exhausted in about 
50 days, while for SN 1999em this happens later by around 50 
days. Accordingly, the luminosity at the CRW phase is main- 
tained by the stored energy during only 22 days for SN 1987A 
and about 75 days for SN 1999em. 

The subsequent run of the bolometric light curves in the 
optimal models of these SNe is entirely determined by the en- 
ergy deposition of gamma rays from the radioactive decays and 
shows the key role played by ^^Ni and ^^Co in powering the 
light curve (Fig. 24). The luminosity of SN 1999em exceeds the 
gamma-ray luminosity from the peak to the radioactive tail. On 
the contrary, the luminosity of SN 1987 A drops abruptly well 
below the gamma-ray luminosity during the adiabatic cooling 
phase, grows to the maximum after the CRW phase in about 
65 days exceeding the latter at day 60, and forms a wide dome 
on the light curve, declining to the radioactive tail and subse- 
quently diving under it. The wide dome on the light curve of the 
SN 1987A is physically equivalent to the phases of the radia- 
tive diffusion cooling, the exhaustion of radiation energy, and 
the plateau tail in the evolution of SN 1999em in the sense that 
it contains the basic features of all these phases. 

At the stage of the radioactive tail, there is a striking dif- 
ference between SN 1999em and SN 1987A: the internal gas 
energy of the envelope decreases with time in SN 1999em and 
increases in SN 1987A (Fig. 23); the bolometric luminosity of 
SN 1999em coincides with the corresponding gamma-ray lu- 
minosity, while the bolometric luminosity of SN 1987A runs 
below it (Fig. 24). The gas energy is controlled by the energy 
rate of gamma-ray deposition, the energy losses due to the adi- 
abatic expansion, and the net rate of absorption-emission pro- 



cesses. At a given phase the last is dominated by the emission 
processes. In the envelope of SN 1999em, the adiabatic energy 
losses are much greater than those in the envelope of SN 1987A 
where the thermal gas energy is much less than that of SN 
1999em. This fact explains both the decrease in the total gas 
energy with time in SN 1999em and its increase in SN 1987A. 

The radiation energy is ruled mainly by the energy rate of 
gamma-ray deposition and the rate of work done by radiation 
pressure. In the envelope of SN I999em the total radiation en- 
ergy is much less than that in the envelope of SN 1987 A, and, 
consequently, the rate of work done by radiation pressure in SN 
I999em is much less than that in SN I987A as well. It leads 
to a negligible difference between the bolometric and gamma- 
ray light curves in the optimal model of SN 1999em and to 
the bolometric light curve running below the gamma-ray light 
curve in the optimal model of SN 1987 A. Note that the higher a 
rate of work done by radiation pressure in an SN envelope, the 
larger the excess of the gamma-ray luminosity over the bolo- 
metric one at the radioactive tail when the SN envelope still 
remains optically thick for gamma rays. In the optimal model 
of SN 1987A, the excess is nearly 7% (Fig. 24), which corre- 
sponds to an underestimation of the ^^Ni amount of the same 
value in equating the observed bolometric luminosity to the 
gamma-ray one. Such an excess of the gamma-ray luminosity 
over the bolometric one is of great interest and should be taken 
into account in measuring the ^^Ni amount from the observed 
bolometric luminosity at the radioactive tail. 

It is very important that the approximation of homologous 
expansion may be used in the atmosphere models starting from 
nearly day 2.8 (Sect. 5.3) for SN 1999em, a normal SN IIP, and 
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from day 1 (Utrobin 2004) for SN 1987A, a pecuUar SN IIP. 
Because of higher expansion velocities in the ejecta, the effects 
of expansion opacity and limb darkening are more prominent 
in SN 1987A than in SN 1999em. For both SNe IIP, the crucial 
role of the time-dependent approach in atmosphere models at 
the photospheric epoch is evident from a comparison of spec- 
tral lines computed in this approach with those in the steady- 
state approximation. 

9. Discussion 

The aim of this paper is to explore the SN 1999em event, a 
normal SN IIP, by comparing the hydrodynamic models and the 
time-dependent atmosphere models with both the photometric 
and spectroscopic observations. The resulting optimal model 
of SN 1999em succeeds in reproducing the observed features. 
The bolometric light curve and the spectral evolution of the 
Ha Une are consistent with a radius of the pre-SN of 5QQRq, a 
mass of the ejected envelope of 19M0, an explosion energy of 
1.3 X 10^' erg, and a mass of radioactive ^^Ni of O.O36M0, the 
bulk of which is confined to layers ejected with velocities less 
than » 660 km s"' . 

The observational data of an SN are measured with some 
errors that yield the uncertainties in values inferred from them. 
In the case of SN 1999em, we can translate these measurement 
errors of the observed parameters into uncertainties of the phys- 
ical parameters using the approximate formulae (2-4). To as- 
sign an uncertainty to the bolometric luminosity in the middle 
of the plateau, we adopt a relative error of 7% for measuring the 
bolometric luminosity itself (Elmhamdi et al. 2003) and an un- 
certainty of ±1 Mpc for the Cepheid distance to the host galaxy 
(Leonard et al. 2003). The resulting relative error in the bolo- 
metric luminosity in the middle of the plateau is 26%, as well 
as the error in the mass of radioactive ^^Ni. The relative errors 
in the photospheric velocities do not exceed 5%, inferred from 
the Fe 11 4924, 5018, and 5169 A absorption lines in the mid- 
dle of the plateau (Hamuy et al. 2001; Leonard et al. 2002). The 
scatter in the explosion date yielded by the EPM varies within 
±2 days (Hamuy et al. 2001; Leonard et al. 2002; Elmhamdi 
et al. 2003). To estimate the measurement errors in the char- 
acteristic duration of the light curve, we consequently adopted 
an uncertainty of +4 days. All the accepted errors in the ob- 
served values yield the following uncertainties in physical pa- 
rameters: +200Rq in the pre-SN radius, ±1.2Mq in the ejecta 
mass, +0.1 X 10^' erg in the explosion energy, and +O.OO9M0 
in the mass of radioactive ^*Ni. 

Elmhamdi et al. (2003) note a clear flattening in the BVRI 
light curves of SN 1999em after the main plateau at the begin- 
ning of the radioactive tail and call it "the second plateau". The 
"UBVRI" bolometric light curve of SN 1999em consequently 
reflects the second plateau. The observed second plateau coin- 
cides in time with the plateau tail of the bolometric light curve 
for the optimal model (Sect. 5.7), but runs slightly below a line 
of the radioactive tail in contrast to the latter. It seems that the 
real bolometric luminosity of SN 1999em is greater than the 
"UBVRI" luminosity during the second plateau. Elmhamdi et 
al. (2003) also report on the second plateau on the tail for SN 
1991G and SN 1997D, and find that its duration is correlated 



with the amount of ejected ^^Ni, as well as the duration of the 
plateau tail does. Of course, a second plateau nature that is dif- 
ferent from that of the plateau tail cannot be excluded at the 
present time. 

In the same study, Elmhamdi et al. explore the width and 
position of the He I 10 830 A Une, which is sensitive to the 
nonthermal ionization and excitation produced by the radioac- 
tive decays, and show that in the ejecta of SN 1999em the 
bulk of ^^Ni was distributed inside a region with a velocity less 
than 1100 km s~^ in the close hemisphere. Obtained in the op- 
timal model, a weak ^''Ni mixing in velocity space up to x 660 
kms~' agrees well with this ^^Ni distribution. By interpreting 
the temporal evolution of the [O I] 6300, 6364 A doublet profile 
in terms of the dust formation, they conclude that the dust oc- 
cupies a sphere with a velocity of » 800 km s ' . Dust grains are 
formed by heavy elements, so we should estimate the volume 
containing them. The characteristic velocity of the outer edge 
of the metal-rich region is nearly 850 kms"' in the optimal 
model. Moreover, the ^^Ni bubble shell containing heavy ele- 
ments is unstable to the Rayleigh-Taylor instability and, con- 
sequently, should be broken by it (Basko 1994) and mixed, 
thereby producing dense clumps in which dust can form. Thus, 
inside the sphere of a velocity of roughly 850 km s"' in the op- 
timal model of SN 1999em, there are conditions that favor the 
dust formation in the ejecta at late times. In addition, from the 
observed [O 1] 6300, 6364 A doublet luminosity, they infer the 
oxygen mass of ~ 0.3 - 0.4Mo, which is consistent with the 
oxygen mass of O.53M0 estimated for the optimal model. 

Baklanov et al. (2005) carry out the hydrodynamic study of 
SN 1999em in multi-group approximation, compare the calcu- 
lated hydrodynamic model to both the observed UBVRI light 
curves and the photospheric velocity, and advocate the hydro- 
dynamic model for the distance of 7.5 Mpc with the pre-SN 
radius of 450 Rq, the ejecta mass of 15 Mq, and the explo- 
sion energy of 0.7 x 10^' erg. However, the low metallicity of 
Z = 0.004 for a red giant pre-SN and the uniform chemical 
composition with X = 0.7 and Z = 0.004 from the surface to 
its center are questionable. It is worth noting that their model 
for the distance of 12 Mpc is not far from our optimal model in 
the basic parameters: the pre-SN radius of 1000 Rq, the ejecta 
mass of 18 Mq, and the explosion energy of 10^' erg. 

The ejected envelope of 19Mq and a neutron star of 
1.58Mo in the optimal model of SN 1999em correspond to a 
2O.58M0 pre-SN star. Assuming an extremely low mass loss, 
we may now conclude that the lower mass limit of the pro- 
genitor of SN 1999em on the main sequence is nearly 21M0. 
Heger et al. (2000) report the results of stellar evolution calcu- 
lations of non-rotating main-sequence stars in the mass range 
of 10 - 25Mq. It turns out that a main-sequence star of 25 M0 
evolves to the pre-SN star of 18.72M0 with the final helium 
core mass of 7.86M0. A linear extrapolation for the pre-SN 
mass of 2O.58M0 yields the stellar mass of 28.9M0 on the 
main sequence and the final helium core mass of 9.6M0. Thus, 
the pre-SN of SN 1999em most likely developed from a main- 
sequence star in the mass range of 21 - 29M0. 

Generally speaking, there is an opportunity to put con- 
straints on the masses and the evolutionary states of the pro- 
genitors of core coUapse SNe from direct observations. But it 
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is a difficult task to identify the progenitors of SNe, and it is 
impossible to predict the time of explosion of a massive star 
precisely. For these reasons Smartt (2002) suggests estimat- 
ing the initial masses of the progenitors by directly identifying 
SNe progenitors in the archive images of outburst sites taken 
prior to the explosions, by estimating the bolometric luminosity 
limits of progenitors as a function of stellar effective tempera- 
ture, and by comparing them with stellar evolutionary tracks 
on the Hertzsprung-Russell diagram. Applying this method to 
SN 1999em, Leonard et al. (2003) have derived the upper mass 
limit for the progenitor of 20 + 5Mq for the Cepheid distance 
to NGC 1637, and Smartt et al. (2003) obtained the upper limit 
of 15Mq for the distance of 11 Mpc. Evidently, our estimate of 
the progenitor mass is in good agreement with the first upper 
mass limit. 

It should be emphasized that the estimates of the bolometric 
luminosity for progenitors are hardly affected by the circum- 
stellar dust, and a comparison with stellar evolutionary tracks 
is far from being unambiguous as the well-known case of SN 
1987 A demonstrated. There is growing observational evidence 
of the crucial role played by circumstellar dust in evaluating 
the bolometric luminosity of a progenitor star. For example. 
Barlow et al. (2005) show that a large fraction of the 5 mag of 
extinction within the host galaxy NGC 6946 toward SN 2002hh 
might be due to circumstellar gas dust, which might condense 
within a stellar wind from an earlier M supergiant or luminous 
blue variable phase of the evolution of the progenitor star. And 
Massey et al. (2005) find that a significant fraction of red su- 
pergiants in Galactic OB associations and clusters show up to 
several magnitudes of excess visual extinction compared to OB 
stars in the same regions and argue that this is very likely due 
to circumstellar dust around the red supergiants. Note that the 
existence of circumstellar dust before an SN outburst does not 
contradict the lack of influence from dust on the photometric 
properties of an SN. Dust grains are most likely evaporated by 
the intense radiation flash during the shock breakout inside a 
sphere of radius of ~ 0.1 pc (Dwek 1983). 

A good fit of the calculated bolometric light curve for the 
optimal model to that observed for SN 1999em is achieved 
with mutual mixing of the hydrogen-rich and helium-rich mat- 
ter at the helium/hydrogen composition interface. The sharp 
boundary between the hydrogen-rich and helium-rich layers 
(Fig. 16a) at the edge of the helium core — a characteristic of 
the evolutionary models of pre-SNe — changes the light curve 
producing a distinct bump at the end of the plateau (Fig. 16b). 
Note that the same feature was obtained by Chieffi et al. (2003) 
in modelling the theoretical light curves of SNe IIP with the 
evolutionary models of pre-SNe. In addition to a deep hydro- 
gen mixing in velocity space downward to « 700 kms~\ the 
optimal model of SN 1999em is characterized by a weak ^*Ni 
mixing up to » 660 km s ' (Fig. 2 and Table 6). 

The hydrodynamic model for SN 1987A is, in turn, char- 
acterized by a moderate ^^Ni mixing in velocity space up to 
~ 3000 kms ' and a deep hydrogen mixing downward to ~ 
600 kms~^ (Utrobin 2005). Moreover, the observations of SN 
1987 A provide clear evidence of moderate mixing of the bulk 
of radioactive ^^Ni up to velocity ~ 3000 kms ' and deep hy- 
drogen mixing down to ~ 500 km s"^ . For example, the [Ni II] 



6.64 jum profile at day 640 gives a velocity vfwhm = 3100 
km s"' and its modelling results in a maximum velocity of 2600 
km s~' (Colgan et al. 1994). The fact that the Ha profile on day 
498 (PhilUps et al. 1990) is not flat-topped implies that there is 
no large cavity free of hydrogen at the center of the envelope 
and that hydrogen is mixed downward to ~ 500 km s'^. 

Within the framework of the neutrino-driven explosion 
mechanism for core collapse SNe, Kifonidis et al. (2003, 2006) 
carried out two-dimensional simulations of explosion models 
with high-order mode perturbations and those with low-order 
(1 = 2 and / = 1) unstable modes induced by the initial global 
deformation of the shock. The low-mode explosion models ex- 
hibit final iron-group velocities of ~ 3300 km s"^ and deep hy- 
drogen mixing downward to a velocity of ~ 500 km s"' , which 
are consistent with the observations of SN 1987A. The moder- 
ate ^*Ni mixing up to a velocity of ~ 3300 km s~^ is the result 
of the larger initial maximum velocities of metal-rich clumps, 
compared to the high-mode models, which protect the fastest 
clumps from the strong interaction with the reverse shock that 
forms below the helium/hydrogen composition interface. The 
initial global deformation of the shock leads to the growth 
of the Richtmyer-Meshkov instability that in turn results in 
the strong inward mixing of hydrogen at the heUum/hydrogen 
composition interface. 

It is Ukely that the strongly anisotropic explosions are able 
to explain both weak and moderate ^^Ni mixing, depending on 
the pre-SN structure and the character of the interaction be- 
tween metal-rich clumps and the reverse shock. Thus, the weak 
^^Ni mixing and the deep hydrogen mixing, which are the char- 
acteristic features of the optimal model for SN 1999em, cannot 
be excluded within the framework of the neutrino-driven ex- 
plosion mechanism (Kifonidis 2006). These promising capa- 
bihties make our interpretation of the SN 1999em event more 
physically grounded and confirm the suggestion of Chieffi et al. 
(2003) that in reality the inner layers of the expeUed envelope 
are strongly mixed during the explosion. 

Based on the hydrodynamic models of SNe IIP, Litvinova 
& Nadyozhin (1983, 1985) constructed simple approximate 
formulae for evaluating the initial radius, the ejecta mass, and 
the explosion energy — three basic physical parameters — 
from the observed properties of individual SN IIP. Using these 
approximate formulae for SN 1999em and assuming the dis- 
tance to the host galaxy to be 11.08 Mpc, Nadyozhin (2003) 
derived the initial radius of 414 Rq, the ejecta mass of 15.0 
Mq, and the explosion energy of 0.68 x 10^' erg. The differ- 
ence in the ejecta mass and the explosion energy, compared to 
those of the optimal model, results mainly from neglecting the 
influence of radioactive ^^Ni on the light curve and using the 
approximation of equilibrium radiative diffusion in their hy- 
drodynamic models. 

We go one step further to construct the approximate formu- 
lae (2-4), which include the influence of radioactive ^^Ni on the 
light curve and are based on the hydrodynamic models com- 
puted in terms of radiation hydrodynamics in the one-group 
approximation with non-UTE effects in the average opacities 
and the thermal emissivity, with nonthermal ionization, and 
with the contribution of lines to opacity. The resulting approx- 
imate formulae, obtained in space of the physical parameters. 
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are valid only for the vicinity of the optimal model because of 
the sensitivity of the bolometric light curve to the pre-SN struc- 
ture and to its chemical composition. To construct the global, 
approximate relations covering a wider range of the parame- 
ters, we have to use a set of hydrodynamic models computed 
for many real SNe IIP. The optimal model for SN 1999em is 
the first element in this set. 

10. Conclusions 

We have presented a comprehensive study of a normal Type 
IIP SN 1999em. Our main aim in this paper was to show the 
absolute necessity of the simultaneous interpretation of both 
the photometric and spectroscopic observations. To do this, we 
calculated the hydrodynamic models of SN 1999em comple- 
mented by the atmosphere models with the time-dependent ki- 
netics and energy balance. We focused mainly on the optimal 
hydrodynamic model of SN 1999em and subsequent compari- 
son to a pecuUar Type IIP SN 1987A. Our results can be sum- 
marized in the following conclusions. 

- The bolometric light curve of SN 1999em and its spectral 
evolution of the Yia line are consistent with a radius of the 
pre-SN of 500 ± 200^0, a mass of the ejected envelope of 
19.0 ± 1.2M0, an explosion energy of (1.3 ±0.1) x 10^^ erg, 
and a mass of radioactive ^^Ni of 0.036 ±0.009Mo, the bulk 
of which is confined to layers ejected with velocities less 
than « 660 km s"^ The optimal hydrodynamic model with 
the helium core mass in the range of 5.6-8.1 Mq matches 
the observed bolometric light curve. 

- The adequate hydrodynamic and atmosphere models of 
SN 1999em distinguish between the short distance of 7.85 
Mpc, the average value of the EPM distance estimates to the 
host galaxy, and the Cepheid distance of 11.7 Mpc. They 
are inconsistent with the short distance of 7.85 Mpc, which 
should be discarded. 

- It is shown that the hydrogen recombination in the atmo- 
sphere of a normal Type IIP SN 1999em, as well as most 
likely other SNe IIP, at the photospheric epoch is essentially 
a time-dependent phenomenon in accordance with the anal- 
ysis of Utrobin & Chugai (2002, 2005). 

- This study reveals the time development of a normal SN IIP 
with the following stages: a shock breakout, an adiabatic 
cooling phase, a phase of cooling and recombination wave, 
a phase of radiative diffusion cooling, an exhaustion of ra- 
diation energy, a plateau tail, and a radioactive tail. 

- The plateau tail of the theoretical bolometric light curve 
coinciding in time with the second plateau observed in SN 
1999em (Ehnhamdi et al. 2003) runs sUghtly above a Une of 
the radioactive tail in contrast to the latter. It seems that the 
real bolometric luminosity of SN 1999em is greater than 
the "UBVRI" luminosity during the second plateau. The 
duration of the plateau tail, as well as the second plateau, is 
correlated with the amount of ejected ^^Ni. 

- The bolometric light curve depends weakly on the chemi- 
cal composition in the outer layers of the ejecta beyond the 
helium core, while hydrogen is abundant there and controls 
the opacity of matter. A mutual mixing of hydrogen-rich 



and heUum-rich matter in the inner layers of the ejecta guar- 
antees a good fit of the calculated bolometric light curve to 
what is observed. Note that the sharp boundary between 
the hydrogen-rich and heUum-rich layers at the edge of the 
helium core — a characteristic of the evolutionary models 
of pre- SNe — changes the hght curve producing an unob- 
served bump at the end of the plateau. 

- A contribution of numerous metal lines to opacity plays a 
fundamental role in reproducing the observed hght curve 
of SN 1999em. The decrease in opacity due to disregarding 
the line contribution speeds up the radiation diffusion and 
causes the bolometric luminosity to increase significantly 
compared to what is observed during nearly the whole out- 
burst and, as a consequence, causes the characteristic dura- 
tion of the light curve to shorten by about 10 days. 

- The comparison of a normal Type IIP SN 1999em with a 
peculiar Type IIP SN 1987 A reveals two very important 
results for SN theory. First, the masses of hehum cores in 
the pre-SNe of SN 1999em and SN 1987A are close enough 
to suppose that in the final stages of stellar evolution nearly 
the same iron cores form within the pre-SNe. This fact and 
roughly the same explosion energies of SN 1999em and 
SN 1987 A together imply a unique explosion mechanism 
for these core collapse SNe. Second, the optimal model for 
SN 1999em is characterized by a weaker ^^Ni mixing in 
velocity space up to ^ 660 kms"' compared to a moderate 
^*'Ni mixing up to ~ 3000 km s"^ in the case of SN 1987A, 
hydrogen being mixed deeply downward to « 700 km s"' 
and ~ 600 kms 

- A significant excess of the gamma-ray luminosity over the 
bolometric one, shown by the optimal model of SN 1987A 
at the radioactive tail, is of great interest and should be 
taken into account in measuring the ^^Ni amount from the 
observed bolometric luminosity at this stage when the SN 
envelope still remains optically thick for gamma-rays. The 
higher the rate of work done by radiation pressure in an SN 
envelope, the larger an excess of the gamma-ray luminosity 
over the bolometric one at the radioactive tail. 

- It is shown that the approximation of homologous expan- 
sion may be used in the EPM and SEAM for determining 
the distances to normal SNe IIP starting from nearly day 3 
after the SN explosion. 

- Based on the hydrodynamic models in the vicinity of the 
optimal model of SN 1999em, we derive the approximate 
relationships between the basic physical and observed pa- 
rameters. The basic physical parameters are the pre-SN 
radius, the ejecta mass, the explosion energy, and the to- 
tal ^^Ni mass. The observed properties of SNe IIP are de- 
scribed by the three parameters: the characteristic duration 
of the light curve, the bolometric luminosity and the pho- 
tospheric velocity in the middle of the plateau. To cover a 
wider range of the parameters, beyond the vicinity of the 
optimal model, we have to construct a global approxima- 
tion that should be based on a set of hydrodynamic models 
of many real SNe IIP. 

Finally, we emphasize that the simultaneous analysis of the 
photometric and spectroscopic observations is in fact crucial 
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for correctly interpreting core collapse SNe and express the 
hope that the obtained basic parameters and properties of SN 
1999em determine its evolutionary state and will be a reliable 
reference for the proper mechanism of SN explosions. 
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